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Abstract 

A novel mechanism to build up a supermassive black hole (SMBH) is proposed. 
' Here, relativistic radiation-hydrodynamic processes are considered from a galactic 

scale to a horizon scale. In the present scenario, the mass of a SMBH is predicted 
to be in proportion to the bulge mass, and it turns out that the SMBH-to-bulge 
mass ratio is basically determined by the nuclear energy conversion efficiency from 
£T) [ hydrogen to helium, e = 0.007. The BH hierarchy is also addressed in relation to the 

BH formation in globular clusters. Futhermore, a new picture for QSO formation 
^ <— | is proposed based on the present SMBH formation scenario. It is predicted that 

a host luminosity-dominant "proto-QSO phase" exists before the QSO phase, and 
q , the proto-QSO phase is preceded by an optically-thick ultraluminous infrared galaxy 

; (ULIRG) phase. 

1.1 Introduction 

The recent discovery of high redshift quasars with z > 6 (Fan et al. 2001) 
implies that the formation of supermassive black holes proceeded in less than 10 9 yr. 

$_i ' Also, the recent compilation of the kinematical data of galactic centers in both active 

' 



and inactive galaxies has shown that a central 'massive dark object' (MDO), which 
is the nomenclature for a black hole (BH) candidate, correlates with the properties 
of galactic bulges. The demography of MDOs have revealed the following relations: 

1) The BH mass exhibits a linear relation to the bulge mass with the ratio of 

/bh = -jT 55 - = 0.001 - 0.006 (1.1) 

J "bulge 

as a median value (Kormendy & Richstone 1995; Richstone et al. 1998; Magorrian 
et al. 1998; Gebhardt et al. 2000a; Ferrarese & Merritt 2000; Merritt & Ferrarese 
2001a; McLure & Dunlop 2002). 

2) The BH mass correlates with the velocity dispersion of bulge stars with a power- 
law relation as Mbh <x cr™, n — 3.75 (Gebhardt et al. 2000b), 4.72 (Ferrarese 
& Merritt 2000; Merritt & Ferrarese 2001a, b), or 4.02 ± 0.32 (Tremaine et al. 
2002). 
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3) The /bh tends to grow with the age of youngest stars in a bulge until 10 9 yr 
(Merrifield et al. 2000). 

4) In disk galaxies, the mass ratio is significantly smaller than 0.01 if the disk stars 
are included (Salucci 2000; Sarzi ct al. 2001). 

5) For quasars the /bh is of a similar level to that for elliptical galaxies (Laor 1998; 
Shields et al. 2002). 

6) The /bh in Seyfert 1 galaxies is not well converged, which may be considerably 
smaller than 0.001 (Wandel 1999; Gebhardt et al. 2000a) or similar to that for 
ellipticals (McLure & Dunlop 2002), while the BH mass-to- velocity dispersion 
relation in Seyfert 1 galaxies seems to hold in a similar way to elliptical galaxies 
(Gebhardt et al. 2000a; Nelson 2000). 

On the other hand, the X-ray emission (Brandt et al. 1997) or Paa lines (Veilleux, 
Sanders, & Kim 1999) intrinsic for active nuclei have been detected in roughly one 
forth of ultraluminous infrared galaxies (ULIRGs). Furhtermore, it has been revealed 
that QSO host galaxies are mostly luminous and well evolved early-type galaxies 
(McLeod & Rieke 1995; Bahcall ct al. 1997; Hooper, Impcy, & Foltz 1997; McLeod, 
Ricke, & Storrie-Lombardi 1999; Brothcrton ct al. 1999; Kirhakos et al. 1999; 
McLure et al. 1999; McLure, Dunlop, & Kukula 2000). All these findings on QSO 
hosts and supermassive BHs imply that the formation of a QSO, a bulge, and a 
supermassive BH is mutually related. 

Some theoretical models have been hitherto considered to explain the BH-to-bulge 
correlations (Silk & Rees 1998; Adams, Graff, & Richstone 2001; Ostriker 2000). But, 
little has been elucidated regarding the physics on the angular momentum transfer 
which is requisite for BH formation, since the rotation barrier by the tidal spin up 
in a growing density fluctuation is given by 



in units of the Schwarzshild radius i?sch, where Mb is the baryonic mass, z is the 
cosmological redshift, and A is the spin parameter which provides the ratio of cir- 
cular velocity to velocity dispersion of dark matter (Umemura 2001). Furthermore, 
required mechanisms for BH formation must work effectively in a spheroidal system 
like a bulge. The a-viscosity or non-axisymmetric gravitational instabilities would 
effectively transfer angular momentum once a disk-like system forms, but they are 
not likely to work in a spheroidal system. 

In this paper, as a potential mechanism that works in a spheroidal system, the 
relativistic drag force by the radiation from bulge stars is considered. As a result, 
the BH-to-bulge ratio is basically determined by the energy conversion efficiency for 
nuclear fusion of hydrogen to helium, e = 0.007. Also, in relation to BH growth, a 
new scenario for QSO formation is addressed. 

1.2 Formation of Supermassive Black Holes 

1.2.1 Mass Accretion due to Radiation Drag 

A radiation hydrodynamic model which could account for the putative corre- 
lations between SMBHs and bulges is recently proposed by Umemura (2001), where 




(1.2) 
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the relativistic drag force by the radiation from bulge stars is considered. The ra- 
diation drag can extract angular momentum from gas and allow the gas to accrete 
onto the center (Umemura, Fukue, & Mineshige 1997, 1998; Fukue, Umemura, & Mi- 
neshige 1997). For the total luminosity L* of a uniform bulge, the radiation energy 
density is given by 

E~L*/cR 2 , (1.3) 

where c is the light speed and R is the radius of the bulge. Then, the angular 
momentum loss rate by the radiation drag is given by 

dlnJ/dt~- X E/c, (1.4) 

where J is the total angular momentum of gaseous component and \ is the mass 
extinction coefficient which is given by \ = with dust absorption coefficient Kd 
and gas density p. Therefore, in an optically-thin regime, 

where r is the total optical depth of the system and M g is the total mass of gas. 
In an optically-thick regime, the radiation drag efficiency is saturated due to the 
conservation of the photon number (Tsuribe & Umemura 1997). Thus, an expression 
of the angular momentum loss rate which includes both regimes is given by 

dlnJ/dt~-^(l-e- r ). (1.6) 

Then, the mass accretion rate is estimated to be 
din J L 
dt 6 

In an optically-thick regime, this gives simply 

M=^, (1.8) 

which is numerically M = 0.1A / /Qyr _1 (i*/10 12 L Q ). This rate is comparable to the 
Eddington mass accretion rate for a black hole with 10 8 M Q , that is, 

M Edd = O.2M yr- V 1 (^gr) . (1-9) 

where rj is the energy conversion efficiency. Unless otherwise stated, rj = 0.42 for an 
extreme Kerr black hole is assumed. The timescale of radiation drag-induced mass 
accretion is 

**- * c i£ - 8 - 6 x ioV << (i<fe) 1 (£) 1 • (lio) 

where i?k P c = -R/kpc and Z is the metallicity of gas. It is noted that the gas which 
is more abundant in metals accretes in a shorter timescale, because the extinction is 
predominantly given by the dust opacity. For the moment, an optically-thick stage is 
considered. Due to the mass accretion induced by the radiation drag, a massive dark 



M = -M fl — ~^(l-e" T ). (1.7) 
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object (MDO) forms at the center of bulge. Then, the mass of MDO is estimated in 
terms of 

A/mdo = / Mdt ~ f ^dt. (1.11) 
Jo Jo c 

In practice, it is likely that optically-thin surface layers are stripped from optically- 
thick clumpy clouds by the radiation drag, and the stripped gas losing angular mo- 
mentum accretes onto the center (Kawakatu & Umemura 2002). 

Next, we employ a simplest analytic model for bulge evolution. The star formation 
rate is assumed to be a Schmidt law, S(t) = kf g . If we invoke the instantaneous 
recycling approximation, the star formation rate is given by 

MjM b = ke~ akt , (1.12) 

where a is the net efficiency of the conversion into stars after subtracting the mass 
loss. The radiation energy emitted by a main sequence star is 0.14e to the rest 
mass energy of the star, m*c 2 , where e is the energy conversion efficiency of nuclear 
fusion from hydrogen to helium, which is 0.007. Thus, the luminosity of the bulge is 
estimated to be 

L % = O.Ueke- akt M h c 2 . (1.13) 

By substituting this in Ql.ll[) . 

Af MDO = O-Uea^M^l - e~ akt ). (1.14) 

The term M(,(l — e~ akt ) represents just the stellar mass in the system which is Mb u i g e 
observationally. As a consequence, the MDO mass to bulge mass ratio is given by 

^ = Q.Wea- 1 = 0.002a i, (1.15) 

-Mbuige 

where ao.5 = a/0.5. It should be noted that the final mass is basically determined by 
e. Also, the proportionality of SMBH mass to bulge mass is a natural consequence 
in the present mechanism. 

1.2.2 Towards the Horizon 

In the above, the prediction for the MDO-to-bulge mass relation is made in 
the context of radiation drag-induced mass accretion. However, the MDO itself is 
not a supermassive black hole, because the radiation drag is not likely to remove the 
angular momentum thoroughly, but some residual angular momentum will terminate 
the radial contraction. Hence, we should consider the further collapse of the MDO 
through other physical mechanisms. 

In the MDO, the viscosity is expected to work effectively because the timescale 
for viscous accretion, j/a v c 2 , is shortened by the radiation drag, where j is the 
specific angular momentum, a v is the viscous parameter, and c s is the sound velocity. 
Thus, the MDO is likely to be a massive self-gravitating viscous disk. For a self- 
gravitating viscous disk, some self-similar solutions are known to give an inside-out 
disk collapse (Mineshige & Umemura 1996; Mineshige & Umemura 1997; Tsuribe 
1999). In particular, Tsuribe (1999) provided a convenient formula for the inside- 
out mass accretion rate, 
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Fig. 1.1. A picture for the formation of a supermassive black hole from a 
galactic scale to a horizon. The radiation from a starburst in a bulge exerts 
the radiation drag onto dusty interstellar gas. Resultantly, the interstellar 
gas sheds angular momentum and accrete onto the center to form a massive 
dark object (MDO). The MDO is equivalent to a massive self-gravitating 
viscous disk. This disk undergoes the inside-out viscous collapse to form a 
rigidly rotating supermassive star with accretion envelope. A rigidly rotat- 
ing supermassive star becomes subject to the general relativistic instability. 
Then, the supermassive star collapse dynamically and eventually a Kerr 
black hole is born. 



M« = (1.16) 

where c s is the sound speed and Q is the Toomre's Q which is kc s /-kGYi for the 
epicycle frequency k and the surface density E. Tsuribe (1999) has found a solution 
of stable accretion with Q « 2. The critical accretion rate is given by 

M a ~ | = O.24M yr- 1 (j^^ ■ (1-17) 

This rate is again comparable to the Eddington mass accretion rate for a black hole 
with 10 8 M Q if T ps 10 4 K [see l|1.9f) ]. Through this inside-out collapse, a central core 
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grows. The core is expected to be a rigidly rotating supermassive star because the 
viscous transfer of angular momentum works to smear out any differential rotation 
in a self-gravitating system. 

The equilibrium configuration and the stability of a rigidly rotating supermassive 
star has been scrutinized by Baumgarte & Shapiro (1999). They found that a rotating 
supermassive star becomes unstable for R < 640/GM/c 2 . As for the dynamical col- 
lapse of a rotating supermassive star, Saijo et al. (2002) performed post-Newtonian 
calculations and found that if a rotating supermassive star is in rigid rotation, it 
can collapse towards the horizon scale without undergoing bar-mode instability. The 
final stage of the collapse was investigated by Shibata & Shapiro (2002) with a full 
general relativistic approach and the emergence of an apparent horizon of a Kerr 
black hole was shown. The resultant spin parameter of the black hole is around 0.75. 

Hence, an MDO formed by radiation drag-induced mass accretion could evolve 
into a supermassive Kerr black hole through the inside-out viscous collapse and the 
general relativistic instability. The present picture for the formation of a supermas- 
sive black hole is summarized in Figure ITTI 

1.2.3 Feedback by AGN Activity 

If the BH accretion causes the nuclear activity, the further radiative mass 
accretion can be induced by the nuclear luminosity Lagn- This feedback works 
positively to grow the SMBH. Consequently, it is predicted that the BH mass to 
bulge mass ratio becomes 

MsH/Mbuige = 0.14£a _1 (l - t?)" 1 = 0.0030^. (1.18) 

This is just comparable to the observed ratio. 

1.2.4 M - a Relation 

If the BH mass is determined by the present mechanism, the BH mass to ve- 
locity dispersion relation is naturally understood in the context of a cold dark matter 
(CDM) cosmology. Supposing the bulge is a virialized system, then GM tot /R = a 2 
and 

R « 0.5i? max oc Af^ 3 (l + w) _1 , (1.19) 

where i? max is the radius at the maximum expansion epoch z max . If a CDM cosmol- 
ogy is assumed, then 

(1 + W) oc M^, (1-20) 

where (3 ~ 1/6 around Mf, = 10 12 M©, almost regardless of the cosmological param- 
eters (Bunn & White 1997). Combining all these relations, we find Mbh oc c n with 
n = 6/(2 — 3/3), which is 4 for j3 = 1/6. This result is just corresponding to the 
inferred relation between the BH mass and the stellar velocity dispersion. 

1.2.5 SMBH in Disks 

The radiation drag efficiency could be strongly subject to the effect of geo- 
metrical dilution (Umemura, Fukue, & Mineshige 1998). If the system is spherical, 
the emitted photons are effectively consumed within the system, whereas a large 
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fraction of photons can escape from a disk-like system and also the drag efficiency 
is considerably reduced across the optically-thick disk. This may be the reason why 
/bh is observed to be significantly smaller than 0.001 in disk systems. 

1.2.6 Hierarchy of Massive BHs 

In the present model, the theoretical upper limit of Jbh is e — 0.007. The 
low mass end of the BH hierarchy is suggested by the recent observations on massive 
BHs in globular clusters, e.g., M B h « 2 x 1O 4 M in Gl (Gebhardt, Rich, & Ho 2002) 
and A/ BH - 3x 10 3 A/ Q in M15 (van der Marel et al. 2002; Gessen et al. 2002). It is 
found that M — a relation still holds in such small spheroidal systems. The present 
radiation hydrodynamic mechanism is expected to work in forming a BH even in a 
small spheroidal system, if the system retains interstellar medium (ISM) sufficient 
to become optically thick. Although in globular clusters a few supernovae can easily 
expel the ISM in an early evolutionary phase, the mass loss from low mass stars in 
later phases may accumulate a sufficient ISM. On the other hand, M33 appears not 
to have a massive BH with an upper limit of 15OOM (Gebhardt et al. 2001). An 
intriguing difference is that M33 contains a significant population of stars younger 
than a few Gyrs (e.g., O'Connell 1983), in contrast to other globular clusters which 
possess massive BHs. Thus, M33 might not have been able to retain a sufficient ISM 
owing to successive supernova explosions. In other words, the BH mass in globular 
clusters may sensitively reflect its star formation history. 

As for the high mass end of the BH hierarchy, whether hypermassive BHs in clus- 
ters of galaxies exist is a crucial test. In galaxy clusters, it is not likely that optically- 
thick intracluster medium is supplied by galactic winds or stripping of ISM. Hence, 
in the present scenario, hypermassive BHs will not form in the centers of galaxy 
clusters. Therefore, the BH hierarchy is predicted to be cut off at the maximum size 
in galaxies, M B h ~ 10 10 M Q . 

1.3 QSO Formation 

Here, we construct a picture of QSO formation based on the present scenario 
for SMBH formation. If the mass accretion driven by the viscosity on to the BH 
horizon is determined by an order of Eddington rate, the BH mass grows according 
to 

M BH = M e" t/tEdd , (1.21) 

where v is the ratio of BH accretion rate to the Eddington rate and iEdd is the 
Eddington timescale, iEdd = 1.9 X 10 8 yr. Mq is the mass of a seed BH with ~ 10 5 M Q 
(Shibata & Shapiro 2002). In order to incorporate the chemical evolution of host 
galaxy, we use an evolutionary spectral synthesis code 'PEGASE' (Fioc & Rocca- 
Volmerange 1997), and also employ a galactic wind model with the wind epoch of 
t w — 7 x 10 8 yr to match the present-day color-magnitude relation (Arimoto & Yoshii 
1987). The system is assumed to change from optically-thick to optically-thin phase 
at t w . Based on the present coevolution model, the evolution of bulge luminosity 
(Lbuige), the AGN luminosity (Lagn), and the mass of SMBH (Mbh) are shown 
in Figure ll~2*l assuming the constant Eddington ratio (y = 1). The Mbh reaches 
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Fig. 1.2. A scenario for QSO formation. The abscissa is time and the 
ordinate is luminosity or mass. L* and Lagn are the bulge luminosity and 
the black hole accretion luminosity, respectively. Mbh is the mass of SMBH. 
The ordinary QSO phase is preceded by a host luminosity-dominant "proto- 
QSO phase" . 



Mmdo at a time i cr0 ss, so that the BH fraction becomes /bh — 0.001, which is just 
comparable to the observed ratio. 

In the optically-thin phase after the galactic wind (t > t w ), Mbh still continues to 
grow until i cross and therefore the AGN brightens with time. After Lagn exhibits 
a peak at i croS s, it fades out abruptly due to exhausting the fuel of the MDO. This 
fading nucleus could be a low luminosity AGN (LLAGN). It is found that the era 
of t w < t < icross can be divided into two phases with a transition time t a it when 
-kbuigc = £agn; the earlier phase is the host luminosity-dominant phase and the 
later phase is the AGN luminosity-dominant phase. The lifetimes of both phases are 
comparable to each other, which is about 10 8 yr. The AGN-dominant phase is likely 
to correspond to ordinary QSOs, while the host-dominant phase is obviously different 
from observed QSOs so far. We define this phase as "a proto-QSO" (Kawakatu & 
Umemura 2003). The observable properties of proto-QSOs are predicted as follows: 
(1) The width of broad emission line is narrower, which is less than 1500km/s. (2) 
/bh rapidly increases from 10~ 53 to 10~ 39 in w 10 8 years. (3) The colors of (B - V) 
at rest bands and ( V — K) at observed bands are about 0.5 magnitude bluer than 
those of QSOs. (4) In both proto-QSO and QSO phases, the metallicity of gas in 
galactic nuclei is Zblr — 8Z Q , and that of stars weighted by the host luminosity is 
~ 3Z , which are consistent with the observations for QSOs and the elliptical 
galaxies. (5) A massive dusty disk (> 1O 8 M ) surrounds a massive BH, and it may 
obscure the nucleus in the edge-on view to form a type 2 nucleus. The predicted 
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Fig. 1.3. Schematic sketch for the coevolution of a SMBH and host galaxy. 

properties of proto-QSOs are quite similar to those of radio galaxies at high redshifts. 
Thus, high-z radio galaxies are a key candidate for proto-QSOs. 

The proto-QSO phase is preceded by a bright and optically thick phase, which may 
correspond to a ultraluminous infrared galaxy (ULIRG) phase. Also, the precursor 
of ULIRGs is an optically-thin and very luminous phase with the lifetime of ~ 10 7 
years. This may correspond to the assembly phase of LBGs or Lya emitters. In this 
phase, the metallicity is subsolar (Z* < O.IZq), and the hard X-ray luminosity is 
L x ~ 5 x 1O 8 L0 if L x = O.ILagn- Such a coevolution scenario of a SMBH and the 
host is summarized in Figure IT31 

1.4 Conclusions 

The SMBH could form through the relativistic radiation-hydrodynamic pro- 
cesses from galactic scale to horizon scale. The resultant mass of a SMBH is in 
proportion to the bulge mass, and the SMBH mass fraction is determined by the 
nuclear energy conversion efficiency from hydrogen to helium, e — 0.007, which is 
the theoretical upper limit in the present mechanism. Massive BHs can form also in 
globular clusters, but hypermassive BHs are unlikely to form in the centers of galaxy 
clusters. As for QSO formation, it is predicted that a host luminosity-dominant 
"proto-QSO phase" exists before the QSO phase. The proto-QSO is a phase of 
growing BH with broad emission line width of less than 1500km/s. The proto- 
QSO phase is preceded by an optically-thick ultraluminous infrared galaxy (ULIRG) 
phase. Before the ULIRG phase, an optically thin starburst phase should exist. 
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